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Abstract 

We present preliminary results on three SLD analyses: the gluon energy spectrum in 
3-jet bbg events, the rate of g — > bb, and the b fragmentation function in Z° decays. The 
gluon energy spectrum, measured over the full kinematic range, is compared with pertur- 
bative QCD predictions. We set new 95% C.L. limits on the anomalous chromomagnetic 
coupling of the b quark: —0.09 < k < 0.06. g h i is measured to be (3.07 ± 0.71 (stat) ± 
0.66 (syst)) xlO -3 . The inclusive B hadron energy distribution in Z° decays is measured 
for the first time over the full kinematic range, using a novel B hadron energy recon- 
struction technique. Several models of b fragmentation including JETSET + Peterson 
are excluded by the data. The average scaled B hadron energy of the weakly decaying 
B hadron is measured to be x B = 0.713 ± 0.005 (stat) ± 0.007 (syst) ± 0.002 (model). 
All three measurements take advantage of the small and stable SLC interaction point as 
well as the excellent vertexing and tracking capabilities of the upgraded CCD-pixel vertex 
detector. 
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1. Introduction 



One of the powerful aspects of QCD is that it provides a framework to understand the 
production of hadronic jets. Perturbative QCD (pQCD) in principle allows us to calculate 
infrared finite quantities as an expansion of the strong coupling a s and its predictions 
are subject to precision experimental tests. The nonperturbative fragmentation process, 
which transforms quarks and gluons into hadrons, however, has not been well-understood 
theoretically, nor has it been well-measured experimentally. This limits the precision in 
our understanding of the production and structure of hadronic jets and will likely affect 
the precision of physics studies at future colliders, where large numbers of hadrons will 
be produced. In view of this, it is important both to test pQCD and to precisely measure 
quantities that will help extract the nonperturbative effects in the fragmentation process. 
Heavy flavors are especially well-suited for this purpose. The mass of the b quark provides 
a natural cutoff in pQCD calculations, a feature that allows a number of precise pQCD 
predictions to be made. Furthermore, e + e _ — > Z° environment is a particularly clean. 
When a B hadron is highly boosted, its decay vertex is significantly displaced from the 
interaction point (IP), allowing high efficiency and purity B tagging. 

Here we report three SLD [|1|] analyses involving heavy quark production. SLD's 
strength to study B physics largely resides in the upgraded CCD-pixel vertex detector 
(VXD3) D. Using VXD3, the impact parameter resolution is 11/im (23/xm) transverse 
to (along) the beam axis for high momentum tracks. SLC's small and stable IP can be 
measured to a resolution of about 4.4/zm in the plane transverse to the beam axis. This 
allows not only the selection of a high purity B sample, but also the precise measurement 
of kinematic quantities. 

2. The Gluon Energy Spectrum in Z° — > bbg 

In studying e + e~ — > 3-jet events , the key is to know the origin of each jet. In previous 
studies, jets were usually energy-ordered and the lowest is assumed to be the gluon jet. 
This is correct about 80% of the time. However, the bias towards selecting low energy 
gluon jets makes this method undesirable for testing pQCD predictions of the gluon jet 
energy spectrum, to which high energy gluon jets are most sensitive. 

In this analysis |3[], the gluon jet is tagged explicitly by identifying the two jets that 
contain B hadrons. A 3-jet event (selected using the JADE algorithm at y cut = 0.02) 
is tagged as a bbg event if exactly two of the three jets each contain two or more tracks 
with normalized transverse impact parameter djo^ > 3. The remaining jet, which does 
not satisfy this criterion, is tagged as the gluon jet. In 2.5% (12.5%) of the selected 
bbg events, the (second) highest energy jet is tagged as the gluon jet, covering the full 
kinematic range. Non-66 events (4%), non-bbg events (21%), and bbg events in which the 
gluon jet is wrongly tagged (1%) are treated as backgrounds and are subtracted from the 
data. The resulting distribution of the scaled gluon energy x g = 2E 9 /Ecm is corrected 
for effects of efficiency and resolution. The fully corrected spectrum is shown in Figure 
1. The cutoff at low x g is caused by the finite y cut value used in jetfinding. Leading order 
(LO) and next-to-leading order (NLO) QCD predictions Q describe the general but not 
the detailed behavior of the spectrum. Higher order QCD effects are clearly important. 
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Parton shower Monte Carlo is able to describe the data. 

The gluon energy spectrum is particu- 
larly sensitive to the presence of an anoma- 
lous chromomagnetic coupling of the b quark 
in the QCD Lagrangian ||, 
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Figure 1. The measured scaled gluon en- 
ergy spectrum. 
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where k and k parameterize the anoma- 
lous chromomagnetic and chromoelectric mo- 
ments, respectively, which might arise from 
physics beyond the SM. Setting k to zero, 
a fit of the theoretical prediction to the data 
yields k = —0.02 ± 0.04 (preliminary), which 
is consistent with zero. The 95% C.L. limits 
are —0.09 < k < 0.06 (preliminary). 



3. The Rate of Secondary bb Production via g 



bb 



Max(Vertex Mass(Jetl), Vertex Mass(Jet2)) 

to [ 

Preliminary 



The process of gluon splitting into a quark- ant iquark pair is poorly known, both theo- 
retically and experimentally, despite the fact that this is one of the elementary processes 
in QCD. The rate of g — > bb, g bb , defined as the fraction of hadronic events in which a 
gluon splits into a bb pair, is an infrared finite quantity and can be calculated using pQCD 
because the mass of the b quark provides a natural cutoff . However, g b i is sensitive to 
both the A|j^ parameter and the b quark mass, which results in a substantial theoretical 
uncertainty in the calculation of g bb . The limited accuracy of the g — > bb prediction is one 

of the main sources of uncertainty in the mea- 
surement of the partial decay width of the 
Z° — > bb. In addition, about 50% of the 
B hadrons are produced via gluon splitting 
at Tevatron, and an even larger fraction is 
expected at the LHC. A better knowledge 
of this process will improve the theoretical 
modeling and predictions of heavy flavor pro- 
duction at such colliders. 

Since the background in this analysis is 
very large, the main task is to enhance the 
signal to background ratio. The 4-parton fi- 
nal state of the signal, qqbb, suggests that 
we select 4-jet hadronic events. We used the 
Durham algorithm with y cut = 0.008. Most 
Z° — > qq events and Z° — > qqg events in 
which the gluon does not split are rejected, 

but the majority of the selected events are still backgrounds from qqgg and qqcc. Back- 
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Figure 2. The larger of the two vertex 



masses. 
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ground of the type qqgg where q 7^ b or c can be mostly removed by looking for heavy 
hadrons in the selected events. Using a topological vertexing algorithm |7j], secondary 
vertices are searched for in the two jets forming the smallest angle. About 300 events in 
which both jets contain a secondary vertex are selected. Figure 2 shows the distribution 
of the larger of the two vertex masses, together with the breakdown of Monte Carlo events 
into the signal and the backgrounds. 

We then focus on distinguishing the topology of g — > bb events from that of bbgg events. 
In many events, the two selected jets actually originate from the same b jet. Figure 3 
shows the angular separation between the two vertices. The high IP and vertex resolutions 
give us a good discriminating power between the signal and the background, about half 

of which peak at cosd ~ 1. In order to en- 
hance the signal to background ratio, we re- 
quire —0.2 < cos6> < 0.96. Since the radi- 
ated virtual gluon in Z° — > qqg is polarized in 
the plane of the three partons, the subsequent 
gluon splitting favors g — > qq emission out of 
this plane. We therefore consider a variable 
similar to the Bengtsson-Zerwas angle, which 
is the angle between the plane formed by the 
two jets containing secondary vertices and the 
plane formed by the other two jets. Figure 4 
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Figure 3. Distribution of the cosine of 
the angle between the two jets forming the 
smallest angle, for data (points) and in 
Monte Carlo (histogram). 
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Figure 4. Distribution of the cosine of 
the angle between the two planes, for data 
(points) and in Monte Carlo (histogram). 



shows the distribution of | cosa |. Events with 
| cosct |> 0.8 are rejected. Finally, since b jets 
from the gluon splitting tend to be softer than 
other jets, we require the total energy of the 
two selected jets to be less than 36 GeV. Af- 
ter these cuts, the background consists of 49% 
bbgg, 34% g — > cc and the remaining 17% 
qqgg, where q 7^ b or c. To further reduce the 
g — > cc background, we require the maximum 
vertex mass to be greater than 2.0 GeV. 62 
events in the data survived all these selection 
cuts. We subtract the expected background of 
27.6±1.2 events and divide the selection effi- 
ciency of 3.86% to obtain the preliminary mea- 
sured value of g b i 



g bl = (3.07 ± 0.71(stat) ± 0.66(syst)) x 10" 



(2) 



4. The b Quark Fragmentation Function 

According to the factorization theorem, the heavy quark fragmentation function can be 
described as a convolution of perturbative and non-perturbative effects. For the b quark, 
the perturbative calculation is in principle understood p, [| . Nonperturbative effects have 
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been parametrized in both model-dependent [LI], [12|, [H| [8] and model-independent 
approaches 0, [15|, [16| . 

It is experimentally challenging to measure the b quark fragmentation function to a 
level of precision sufficient to distinguish among the various models. Since the b quark 
fragmentation function is the probability distribution of the fraction of the momentum 
of the b quark carried by the B hadron, the most sensitive experimental determination 
of the shape of the b fragmentation function is expected to come from a precise direct 
measurement of the B hadron energy (or momentum) distribution. The difficulty of such 
a measurement stems mostly from the fact that most of the B decays can only be partially 
reconstructed, causing a fraction of the B energy to be missing from the B decay vertex. 



Recent direct measurements at LEP [H], |18| and SLD [19| have used overall energy- 



momentum constraints and calorimetric information to extract this missing energy in a 
sample of semi-leptonic B decays. These measurements suffer from low statistics as well 
as poor B energy resolution at low energy, and hence have a relatively weak discriminating 
power between different shapes of the fragmentation function. Indirect measurements [p0[] 
such as the measurement of the lepton spectrum have been used to constrain the average 
B energy. These measurements, however, are not sensitive to the shape of the B hadron 
energy distribution. 

Here we describe a new method for reconstructing individual B hadron energy with 
good resolution over the full kinematic range while achieving a much higher efficiency |2l|] . 
We reconstruct secondary B decay vertices inclusively with high efficiency using a topo- 
logical vertexing algorithm At SLD, the B flight direction, pointing along the line 
joining the IP and the secondary vertex, is well-measured because of the very small IP 
error and the excellent vertex resolution. Therefore the transverse momentum P t of tracks 
associated with the vertex relative to the B flight direction is also well-measured. We then 
obtain the invariant mass M c h and the total energy of the associated tracks, assigning 
the pion mass to each charged track. Furthermore, we assume that the true mass of the B 
hadron decayed at the vertex, Mb, is equal to the B° meson mass. An upper limit on the 
mass of the missing particles is then found to be M$ max = M| - 2M B \jM 2 ch + P? + M c \. 
Since the true missing mass M^ ue is often rather close to M 0max , M 0max is subsequently 
used as an estimate of M^~ ue to solve for the longitudinal momentum of the missing parti- 
cles from kinematics, and hence the energy of missing particles Eq. The B hadron energy 
is then Eb = E c h + Eq. Since < M^ ue < M 0max , the B energy is well-constrained 
when Mo max is small. In addition, most uds and c backgrounds are concentrated at large 
M 0max . We choose an ad hoc upper cut on the M§ max to achieve a nearly xb -independent 
B selection efficiency of about 3.9%; the estimated purity is about 99.5%. A total of 1920 
vertices in the 1996-97 data satisfy all selection cuts. Figure 5 shows the distribution of 
the reconstructed scaled weakly decaying B hadron energy for data and Monte Carlo. 

We examine the normalized difference between the true and reconstructed B hadron 
energies for Monte Carlo events. The distribution is fitted by a double Gaussian, resulting 
in a core width (the width of the narrower Gaussian) of 10.4% and a tail width (the width 
of the wider Gaussian) of 23.6% with a core fraction of 83%. The core width is essentially 
independent of xb, another feature that makes this method unique. 

After background subtraction, the distribution of the reconstructed scaled B hadron 
energy is compared with a set of ad hoc functional forms of the xb distribution in order 
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to estimate the variation in the shape of the xb distribution. For each functional form, 
the default SLD Monte Carlo is re-weighted and then compared with the data bin-by-bin 
and a \ 2 is computed. The minimum \ 2 is found by varying the input parameter (s). 
The Peterson function |T2[ , two ad hoc generalizations of the Peterson function |Tj| , and 
an 8th-order polynomial are consistent with the data. We exclude the functional forms 
proposed by BCFY |L£]], Collins and Spiller f["4jl , Kartvelishvili flOf , Lund [[H| and a power 
function of the form f(x) 
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Figure 5. Distribution of the recon- 
structed scaled B hadron energy. 



Figure 6. The binwise average of the seven 
unfolded distributions (preliminary). 



We then test several heavy quark fragmentation models. Since the fragmentation 
functions are usually functions of an experimentally inaccessible variable (e.g. z = 



(E+p\\)h/ (E+p\\)q), it is necessary to use a Monte Carlo generator such as JETSET [^2 



to generate events according to a given input heavy quark fragmentation function. The 
resulting B energy distribution is then used to re-weigh the Monte Carlo distribution 
before comparing with the data. The minimum x 2 is found by varying the input param- 
eter (s). Within the context of the JETSET Monte Carlo, Kartvelishvili [ITU] , Bowler fTI 



and the Lund [T3| models are consistent with the data, while BCFY [16[, Collins and 



Spiller jTJ], and Peterson [12 1 models are found to be inconsistent with the data. 

The four functional forms and the three fragmentation models consistent with the 
data are then used in turn to calculate a model-dependent unfolding matrix which is 
applied to the data to obtain the unfolded xb distribution. The resulting seven unfolded 
xb distributions show substantial model- variations. Figure 6 shows the binwise average 
of the seven unfolded distributions, where the inner error bar represents the statistical 
error and the outer error bar is the quadrature sum of the r.m.s. of the seven unfolded 
distributions and the statistical error within each bin. The outer error bars therefore 
provide an envelope within which the true xb distribution is expected to vary. The mean 
of the scaled weakly decaying B hadron energy distribution is obtained by taking the 
average of the means of the seven functions. The r.m.s. of the seven means is regarded as 
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the error on model-dependence. The preliminary result is 



< x B >= 0.713 ± 0.005(stat) ± 0.007{syst) ± 0.002(modeZ) (3) 

where the small model-dependence error indicates that < xb > is relatively insensitive to 
the allowed forms of the shape of the fragmentation function. 
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